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1^ ■ We calculate the cross section of J/ijj plus prompt-photon inclusive production 

. in 77 collisions at next-to-leading order within the factorization formalism of nonrel- 

\ ativistic quantum chromodynamics (NRQCD) focusing on direct photoproduction. 

' Apart from direct J/ij: production, we also include the feed-down from directly- 

produced Xcj and tp' mesons. We discuss the analytical calculation, in particular 
Oh! the treatment of the various types of singularities and the NRQCD operator renor- 

malization, in some detail. We present theoretical predictions for the future e^e~ 
linear collider TESLA, taking into account both brems- and beamstrahlung. 
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1 Introduction 



Since the discovery of the J /ip meson in 1974, charmonium has provided a useful labora- 
tory for quantitative tests of quantum chromodynamics (QCD) and, in particular, of the 
interplay of perturbative and nonperturbative phenomena. The factorization formalism 
of nonrelativistic QCD (NRQCD) 1 2J provides a rigorous theoretical framework for the 
description of heavy- quarkonium production and decay. This formalism implies a separa- 
tion of short-distance coefficients, which can be calculated perturbatively as expansions in 
the strong-coupling constant a^, from long-distance matrix elements (MEs), which must 
be extracted from experiment. The relative importance of the latter can be estimated 
by means of velocity scaling rules; i.e., the MEs are predicted to scale with a definite 
power of the heavy-quark [Q) velocity v in the limit f -C 1. In this way, the theoretical 
predictions are organized as double expansions in as and v. A crucial feature of this for- 
malism is that it takes into account the complete structure of the QQ Fock space, which 
is spanned by the states n = ^s+i^^W ^j^j^ definite spin S, orbital angular momentum L, 
total angular momentum J, and color multiplicity a = 1,8. In particular, this formalism 
predicts the existence of color-octet (CO) processes in nature. This means that QQ pairs 
are produced at short distances in CO states and subsequently evolve into physical, color- 
singlet (CS) quarkonia by the nonperturbative emission of soft gluons. In the limit u — > 0, 
the traditional CS model (CSM) |2] is recovered. The greatest triumph of this formalism 
was that it was able to correctly describe j3] the cross section of inclusive charmonium 
hadroproduction measured in pp collisions at the Fermilab Tevatron j3] , which had turned 
out to be more than one order of magnitude in excess of the theoretical prediction based 
on the CSM. Apart from this phenomenological drawback, the CSM also suffers from 
severe conceptual problems indicating that it is incomplete. These include the presence 
of logarithmic infrared (IR) singularities in the 0{as) corrections to P-wave decays to 
light hadrons and in the relativistic corrections to S-wave annihilation [HI, and the lack 
of a general argument for its validity in higher orders of perturbation theory. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and uN deep-inelastic 
scattering (DIS), e~^e~ annihilation in the continuum, Z-boson decays, 77 collisions, and 
6-hadron decays may be found in the literature; for reviews, see Ref. jTj. Furthermore, 
the polarization of ip' mesons produced directly and of J/ip mesons produced promptly, 
i.e., either directly or via the feed-down from heavier charmonia, which also provides a 
sensitive probe of CO processes, was investigated. Until recently, none of these studies 
was able to prove or disprove the NRQCD factorization hypothesis. However, HI data 
of ep ^ e + J/ip + X in DIS at the DESY Hadron Electron Ring Accelerator (HERA) 
[HI and DELPHI data of 77 ^ J/ijj + X at the CERN Large Electron Positron Col- 
lider (LEP2) P provide first independent evidence for it by agreeing with the respective 
NRQCD predictions [TOITT] . 

The verification of the NRQCD factorization hypothesis is presently hampered both 
from the theoretical and experimental sides. On the one hand, the theoretical predic- 
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tions to be compared with existing experimental data are, apart from very few exceptions 
|12|13|14|15j . of lowest order (LO) and thus suffer from considerable uncertainties, mostly 
from the dependences on the renormalization and factorization scales and from the lack of 
information on the nonperturbative MEs. On the other hand, the experimental errors are 
still rather sizeable. The latter will be dramatically reduced with the upgrades of HERA 
(HERA II) and the Tevatron (Run II) and with the advent of CERN LHC and hope- 
fully a future e^e~ linear collider (LC) such as the TeV-Energy Superconducting Linear 
Accelerator (TESLA), which is presently being designed and planned at DESY and has 
just been endorsed by the International Committee for Future Accelerators (ICFA). On 
the theoretical side, it is necessary to calculate the next-to-leading-order (NLO) correc- 
tions to the hard-scattering cross sections and to include the effective operators which are 
suppressed by higher powers in v. 

In this paper, we concentrate on the inclusive production of J/ip mesons in high- 
energy 77 collisions. As mentioned above, this process was studied at LEP2 where 
the photons originated from hard initial-state bremsstrahlung. At high-energy e+e~ LCs, 
an additional source of hard photons is provided by beamstrahlung, the synchrotron 
radiation emitted by one of the colliding bunches in the field of the opposite bunch. 
The highest possible photon energies with large enough luminosity may be achieved by 
converting the e^e~ LC into a 77 collider via back-scattering of high-energetic laser light 
off the electron and positron beams. 

In order for & J/ip meson to acquire finite transverse momentum (pt), it must be 
produced together with another particle or a hadron jet (j). Recently, we studied the 
process 77 J/ip + j + X, where X denotes the hadronic remnant possibly including 
a second jet, at NLO ^3]. In this paper, we perform a similar analysis for the process 
77 — >• J/ip + 'J + X , where 7 represents a prompt photon. This process leads to a 
spectacular signal. On the one hand, J/ip mesons can be easily identified through their 
decays to e~^e~ or ^~ pairs, with a combined branching fraction of B{J/tp —>■ = 
(11.81 ±0.14)%. On the other hand, isolated energetic photons can be detected with high 
efficiency through the distinctive showers they produce in the electromagnetic calorimeter, 
which can be distinguished from those due to 7r° mesons on a statistical basis, or through 
their conversion to e+e~ pairs. 

The incoming photons can interact either directly with the quarks participating in 
the hard-scattering process (direct photoproduction) or via their quark and gluon content 
(resolved photoproduction). Thus, the process 77 —> J/ip + 1 + X receives contributions 
from the direct, single-resolved, and double-resolved channels. All three contributions are 
formally of the same order in the perturbative expansion. This may be understood by 
observing that the parton density functions (PDFs) of the photon have a leading behavior 
proportional to a InlM"^ / Aqqj)) oc a/as, where a is the fine-structure constant, M is the 
factorization scale, and Aqcd is the asymptotic scale parameter of QCD. 

In the following, we focus our attention on the direct channel, for which there is a 
phenomenological and a theoretical argument. The former is exposed already at LO. In 
fact, 77 — > J/ip + 'J + X almost exclusively proceeds through direct photoproduction 
because it is then a CS process jl6|17j . By contrast, the single-resolved contribution 
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is strongly suppressed, by the ratio iO'^l'^ ^Sf^ ) / UD^I"^ 



OC V . There is also 



a CS process in the double-resolved channel. However, it is mediated at the parton 
level via gluon fusion and thus suppressed by the smallness of the gluon PDF of the 
photon. The theoretical argument is related to the observation that 77 — J/V' + 7 + 
via direct photoproduction does not yet require mass factorization at NLO. Therefore, 
this contribution does not suffer from an artificial M dependence, so that it is justified 
to consider it separately. By the same token, the M dependence of the single-resolved 
contribution at LO is relatively feeble and to be compensated by the direct contribution 
at next-to-next-to- leading order (NNLO), which is beyond the scope of our analysis. The 
NLO treatments of the single-resolved and double- resolved channels are left for future 
work. 

In the parton model, prompt photons can also be produced through the fragmenta- 
tion of quarks and gluons, which necessitates the introduction of photon fragmentation 
functions. From arguments similar to those regarding the photon PDFs mentioned above, 
it follows that also the photon fragmentation functions are formally of order a /as- Thus, 
consistency requires that their contribution be taken into account in the theoretical de- 
scription of prompt-photon inclusive production. However, if the prompt-photon momen- 
tum is integrated over, then the enhancement factor ln(M^/AQQj-,) oc l/a^ disappears, so 
that, in direct photoproduction, photon fragmentation is degraded to a NNLO contribu- 
tion. Strictly speaking, we consider the process 77 — J/'?/' + X^, where is a hadronic 
system containing a prompt photon, rather than the process 77 J/?/' + 7 + X. 

In order to distinguish prompt photons from bremsstrahlung photons, some hardness 
condition must be introduced, e.g., by requiring that the transverse momentum of the 
photon be in excess of some minimum value. From the theoretical point of view, this is 
essential because soft-photon bremsstrahlung generates an IR catastrophe, which must be 
removed by including also the corresponding virtual-photon contribution. However, this 
would then be equivalent to calculating the 0{a) correction to the process 77 J /ip + j, 
which is suppressed against the 0{as) one already provided in Ref. [^, and this is not 
our intention here. 

The J/ip mesons can be produced directly; or via radiative or hadronic decays of 
heavier charmonia, such as XcJ and ip' mesons; or via weak decays of h hadrons. The 
respective decay branching fractions are -B(Xco ~^ J/i' + l) = (1-18 ± 0.14)%, B{xci 
J/^ + 7) = (31.6 ± 3.3)%, B{xc2 ^ W + 7) = (20.2 ± 1.7)%, B{i,' ^ J/i, + X) = 
(57.6 ± 2.0)%, and B{B J/ip + X) = (1.094 ± 0.032)% The h hadrons can be 

detected by looking for displaced decay vertices with dedicated vertex detectors, and the 
J/ip mesons originating from their decays can thus be treated separately. Therefore and 
because of the smallness of B{B J/ip + X), J/ip production through 6-hadron decay is 
not considered here. The cross sections of the four residual indirect production channels 
may be approximated by multiplying the direct-production cross sections of the respective 
intermediate charmonia with their decay branching fractions to J/ip mesons. 

To summarize, the goal of the present analysis is to calculate the inclusive cross section 
of 77 J /ip+X-^ in direct photoproduction for finite values of px at NLO within NRQCD 
allowing for the J/ip meson to be promptly produced. The LO result, also including the 
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Table 1: Values of k in the velocity-scaling rule (^0^[n]J oc v'^ for the leading cc Fock 
states n pertinent to H = J/^ip, Xcj, i^'- 



single- and double-resolved contributions, may be found in Refs. |16|17j and the related 
references cited therein. The leading relativistic correction, which originates from the 
cc Fock state n = ^Pi^^ and is of 0{v'^) relative to the LO result, may be extracted 
from Ref. ^Hl- Apart from being of general phenomenological relevance, the analyses 
reported here and in Ref. ^S] should also be of conceptual interest for the theoretical 
heavy- quarkonium community because this is the first time that the full NLO corrections 
are evaluated for inclusive 2^2 processes^ within the NRQCD framework. Several 
theoretical problems that need to be tackled in our present analysis were encountered in 
similar form in Ref. ^H] and do not need to be explained again in detail. However, there are 
also qualitatively new features, especially in connection with the operator renormalization 
and the cancellation of IR singularities, which need to be discussed carefully. 

This paper is organized as follows. In Section |2l we describe our analytical calculation 
in some detail, focussing on the qualitatively new features. In Section El we present our 
numerical results appropriate for the e"'"e~ mode of TESLA, and discuss their phenomeno- 
logical implications. Our conclusions are summarized in Section |3] 



2 Analytic results 



We start this section with a few general remarks. In our analytic calculation, we take 
the color gauge group to be SU(A^c) with a generic value of Nc, which is put equal to 
3 in our numerical analysis. Color factors appearing in our formulas include Tp = 1/2, 
Cf = [N^ — 1) /{2Nc), and Ca = Nc- We work in the fixed-fiavor-number scheme, with 
Hf = 3 active quark flavors q = u,d,s, which we treat as massless. The charm quark c 
and antiquark c, with mass m, only appear in the final state. We denote the fractional 
electric charge of quark g by e^. 

The cc Fock states contributing at LO in v are specified for H = J/ip, Xcj, i^' in TablelH 
Their MEs satisfy the multiplicity relations 



3p{8) 



(2J+ 1) (C"^("5) 



3p(8) 



^Here, the cc bound state is considered as one particle. 
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(2J + 1)(C»^=" ^Sl 



3p(l) 
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which follow to LO in v from heavy-quark spin symmetry. 

We employ dimensional regularization with c? = 4 — 2e space-time dimensions to handle 
the ultraviolet (UV) and IR singularities, and we introduce unphysical 't Hooft mass scales 
IX and A to ensure that the renormalized strong-coupling constant and NRQCD MEs have 
the same mass dimension as their LO counterparts. We formally distinguish between UV 
and IR poles, which we denote as 1/euv and l/ejR, respectively. We apply the projection 
method of Refs. |14|19j . which is equivalent to the dimensional matching procedure 
of Ref. [201 ! order to extract the short-distance coefficients that multiply the MEs. 
However, in order to conform with common standards, we adopt the normalizations of 
the MEs from Ref. |2] rather than from Refs. |13|14j : i.e., the MEs include spin and color 
average factors. 

There is only one partonic subprocess at LO, namely 



7(^1) +7(^2 



cc 



(2) 



where the four-momentum assignments are indicated within the parentheses. The cor- 
responding Feynman diagrams are depicted in Fig. ^ In the NLO analysis, we need to 
evaluate the cross section of process (0) in d space-time dimensions and retain terms of 
0{e) because UV counterterms appear in multiplicative renormalization. The final result 
is listed in Eqs. (1) and (2) of Ref. [^. 

The analogous process with n = ^P^^'' yields a relativistic correction of O(f^), and its 
cross section may be obtained from Eq. (3) of Ref. ^16j by substituting 



2N^CFa, 



O 



H 



lp(8) 



m - 1 



O 



H 



'Pi 



(1) 



(3) 



2.1 Virtual corrections 

The diagrams that generate the virtual corrections to the cross section of process (j2I) are 
obtained by attaching one virtual gluon line in all possible ways to the tree-level seed 
diagrams of Fig. ^ They include the self-energy, triangle, box, and pentagon diagrams 
which are depicted in Figs.|2ta), (b), (c), and (d), respectively. Loop insertions in external 
charm-quark lines are accommodated in the corresponding wave-function renormalization 
constant and are not displayed in Fig. |21 The self-energy and (non-vanishing) triangle 
diagrams are UV divergent; the box diagrams are finite; and the pentagon diagrams are 
IR divergent. The latter also contain Coulomb singularities, which are cancelled after 
taking into account the corresponding corrections to the operator (^O^ ^S['^ V In the 
practical calculation, the Coulomb singularities are first regularized by an infinitesimal 
gluon mass. This regularization prescription is then transformed into one implemented 
with a small relative velocity v between the c and c quarks |12p21j . In contrast to the case 



6 



of 77 ~^ J/-?/' + j + X considered in Ref . , diagrams containing a quark loop vanish 
by color conservation. Notice that, at NLO, virtual corrections only occur in connection 
with the cc Fock state n = 

We apply two independent approaches to calculate the one-loop diagrams. The first 
one uses FeynArts |22] to generate the diagrams and self- written Mathematica codes to 
apply the projectors and provide expressions, which are afterwards treated with a FORM 
program to perform the tensor reduction and the extraction of the UV and IR singularities. 
The result is then transformed into a Fortran code to be used for the numerical evaluation. 
The second approach utilizes QGRAF [23j for the generation of the diagrams, FeynCalc 
for the tensor reduction, and LoopTools [22] for the numerical evaluation of the IR-safe 
integrals. These packages are surrounded by self-written interface programs, which allow 
for a completely automated computation flE]. Our analytic result is too lengthy to be 
presented here. 

The self-energy and triangle diagrams of the type indicated in Fig. |21 contain UV 
singularities, which are removed by renormalizing the charm-quark mass m and wave- 
function ip in the LO diagrams of Fig. ^ We renormalize m and in the on-mass-shell 
scheme. The relevant renormalization transformations may be found in Ref. |15j . 

A crucial feature of effective field theories, such as NRQCD, is that the composite 
operators are generally subject to renormalization. In the case of NRQCD, this is essential 
in order to ensure the complete cancellation of IR and Coulomb singularities at NLO and 
so to overcome the conceptual problems of the CSM mentioned in Section ^ To be 
consistent with the rest of our calculation, we also employ dimensional regularization 
here. We adopt the technique described in Refs. |2|14j to directly evaluate the NLO 
corrections to the NRQCD operators. In this way, we avoid having to match partonic 
cross sections evaluated in NRQCD with their counterparts in full QCD. 

In the case under consideration, we have to renormalize the CS ME (^O^ 
which appears at LO. In d space-time dimensions, this ME has mass dimension d — 
The 't Hooft mass scale of NRQCD, A, is introduced to keep its renormalized version, 
which we wish to extract from experimental data, at mass dimension 3. 

The four-quark operator (^O^ '^Si'^ ^ is related to the amphtude for the elastic scat- 
tering of a cc pair. The corresponding tree- level diagram is depicted in Fig. ^a). The 
one-loop corrections to this amplitude are obtained by attaching a virtual gluon line in all 
possible ways to the external heavy-quark legs, and they involve self-energy and vertex 
corrections [see Fig. |21^b)-(e)]. Using the NRQCD Feynman rules in the quarkonium rest 
frame, expanding the one- loop integrands as Taylor series in 1/m, and performing the 
integration over the loop momentum, we obtain the unrenormalized one-loop result 



H 



1 ^ 2 



XI— --lE^^r^r , (4) 

where the subscript labels the tree-level quantity, /i is the 't Hooft mass scale of QCD 
that enters through the d-dimensional loop integration, and 'Je is Euler's constant. The 
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term proportional io 1/v represents the Coulomb singularity, which arises from the ex- 
change of a longitudinal gluon between the outgoing c and c quarks [see Fig. Efc)]. Obvi- 
ously, NRQCD operators of different cc Fock states n start to mix at one loop. Further- 
more, the presence of UV singularities indicates that they need renormalization. In the 
following, we choose the MS scheme for that. We thus write 



O 



H 



J=0 



(^) + ^^ZZ3. I I exp(-e7i5) 



H 



3p(8) 



(5) 



where the subscript r labels the renormalized quantity and we identify A with the NRQCD 
renormalization scale. Inserting Eq. Q into Eq. (@)), we obtain 



H 



)^ (A) ( 1 



J=0 



3p(8) 



2v 



2iCA-2CF)a, fATTfx'Y 1 
W exp(-e7.)- 



eiR 



(6) 



Substituting Eq. ^ into the LO result generates an IR counterterm at 0{as) that is 
indispensible to render the total NLO result finite. The Coulomb singularity present in 
Eq. © is necessary to cancel similar terms in the virtual corrections. 

The renormalization group equation that determines the A dependence of 



H 



'^S[^^ ) (A) may be derived by differentiating Eq. © with respect to A and then 



taking the physical limit e — > 0. It reads 



A^d 



2{Ca - 2C^)a,(^ ^ /^^ 



3p(8) 



(7) 



There is no obvious physical reason to distinguish between the scales n and A, which refer 
to the ccg and cede vertices in the one-loop diagrams of Figs. infb)-(e), respectively. Both 
scales should essentially be of 0{m). In the following, we thus identify ^ = X. Integration 
of Eq. (|7j) then yields 



H 



(A) 



H 



'Sf' (Ao) + 



S{Ca-2Cf) lj^a.(Ao 



O 



H 



3p(8) 



(Ao) is assumed to be 



where Aq is a reference scale for which the value of (O^ 
known and (3q = (11/3)C^ — (4/3)Tpnj is the one-loop coefficient of the QCD beta 
function. In want of a genuine NLO determination of the MEs (^0^[n]Y (Aq) from a 

fit to experimental data, we choose Xq = m and identify (^0^[n]'^ (m) with their A- 
independent LO values, which are known from the literature |2Zj. Since the LO cross 
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section is of 0{a'^), we have to employ in Eq. (jHJ the one-loop formula for as{X), which 
reads 

a,{X) = ^. (9) 

/5oln(AVA^CD) 

2.2 Real corrections 

The real corrections to the cross section of process Q arise from the partonic subprocess 

'j{h)+-f{h) ^ cc[n]{p)+g{h)+-f{p^), (10) 

where n = ^Sq^\^Pj^\ The corresponding diagrams emerge by attaching one real gluon 
line in all possible ways to the tree-level seed diagrams of Fig. ^ and are presented in 
Fig. m Process (fTUI) with n = '^Si^\^Pj^^ is prohibited by color conservation, since the 
cc pair must be in a CO state to neutralize the color of the gluon. On the other hand, 
process (ITnjl with n = ^5*1^^ is forbidden by Furry's theorem ||2Si, as may be understood 
by observing that the ^S[^^ projector effectively closes the charm-quark line and acts like 
a vector coupling, so that we are dealing with a closed fermion loop containing five vector 
couplings. This was also verified by explicit calculation. 

Generically denoting the transition matrix (T) elements of process fllUp by %., its cross 
sections may be evaluated as 

dar = ^dPS3{h + k2;p,h,p^)\%:^, (11) 
^ s 

where s = {ki + ^2)^ is the square of the center-of-mass energy. Here and in the following, 
we denote the Lorentz-invariant A^-particle phase-space element in d dimensions as 



\ (2^) 



i=l 



(12) 

where pi and rrii are the four-momenta and masses of the outgoing particles and P is the 
total four-momentum of the incoming particles. 

Integrating Eq. (fTT|) over the three-particle phase space while keeping the value of pt 
finite and constraining the final-state photon to be hard, we encounter IR singularities 
of the soft type in the case oi n = ^Pj ■ In fact, in the soft-gluon limit, 0, |7^P 

factorizes as 

CT = Wi^3, (13) 
where Tq is the T-matrix element of process Q and 



(14) 



is the appropriate Eikonal factor. Since the soft gluon does not affect the kinematics of 
the hard process involving the residual particles, the three-particle phase space can be 
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decomposed as 



dFSsih + k2;p,h,p'') = dFS2{h + k2;p,p^)dFSl (15) 

where 

The integration of the Eikonal factor over the soft part of the phase space, can be per- 
formed analytically and yields j29] 

' ' exp(-e7B) - + - In 3 - 77 + C'(e) , (17) 



37rm2 \6^s J \e p 1- 13 3 

where 6 is an infinitesimal dimensionless phase-space-slicing parameter and 
= (s - 4m2) / (s + 4m2). 
Since, in the diagrams of Fig. EJ all photons and gluons are coupled to the massive c- 
quark line, IR singularities of the collinear type do not occur. The integration of Eq. (fTT|) 
over the complementary region of phase space can be performed numerically in d = 4 
space-time dimensions because there are no UV or IR singularities. 



2.3 Assembly of the NLO cross section 

The NLO result for the cross section of process © is obtained by adding the virtual and 
real corrections to the LO result. Collecting the various contributions discussed above, 
arising from the virtual corrections (vi), the parameter and wave- function renormalization 
(ct), the operator redefinition (op), the soft-gluon radiation (so), and the hard-gluon 
emission (ha), we can schematically write the resulting differential cross section of 77 — >■ 
H + X^, including the MEs, as 

da{fi, A) = rfo-o(A)[l + 5vi(/i; euv, em, v) + 4t(yu; euv, em) + Sopifi, A; ejR, v)] 

+ rfaso(/i, A; eiR, S) + da^i^if^, (18) 

where the dependences on the unphysical mass scales /i and A, and the regulators eiR, 
euv, "^^5 and 6 are indicated in parentheses for each term. Everywhere in Eq. (fTHj) . ^^(/i) 
is evaluated in the MS scheme using the one- loop formula Q, is defined in the OS 
scheme, and the MEs are understood as their MS values {0^[n\)r{\). This is necessary 
to ensure the exact cancellation of the singularities. 

The right-hand side of Eq. ()18|) is manifestly finite. The UV divergences cancel between 
(5vi and 6ct', the IR singularities among 6^1, Set, ^op, and daso', and the Coulomb singularities 
between 6^1 and Sop. Note that 6op is UV finite upon operator renormalization; therefore, 
6op in Eq. (fT8|l does not depend any more on euv- The full details about the cancellation of 
the IR and Coulomb singularities in the various cc Fock states n are presented in Table 01 
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Subprocess 



Source 



IR- or Coulomb- singular term 



7 + 7 



7 + 7 



cc 



cc 



3p(8) 



+ 7 



+ ^ + 7 



virtual 
operator 
operator 
soft 



2v \ 



H 



3p(8) 



2(CA-2CF)a, 



Table 2: Compilation of the IR- and Coulomb- singular terms arising from the various 
sources in the various partonic subprocesses. 



The right-hand side of Eq. ()18|) is independent of the cut-off parameter S. The S 
dependence cancels between daso and da^a- While dago is known in analytic form, the 
phase-space integrals occurring in the evaluation of da^s, are rather cumbersome and 
are thus solved numerically. Consequently, the cancellation of 6 has to be established 
numerically. The goodness of this cancellation is assessed in Section El 

The right-hand side of Eq. (fTH|) is independent of the renormalization scales fi and 
A up to terms that are formally beyond NLO; this cancellation is not exact because the 
running of as{fi) and {0^[n])r{X) is determined from the respective renormalization group 
equations, which resum logarithmic corrections to all orders. 



3 Numerical results 

We are now in a position to present our numerical analysis of the inclusive production 
of prompt J/ip mesons in association with prompt photons in two-photon collisions with 
direct photon interactions at NLO in the NRQCD factorization framework. We consider 
TESLA in its e~^e~ mode with y/s = 500 GeV, where the photons are produced via 
bremsstrahlung and beamstrahlung. 

We first specify our input parameters. We use m = 1.5 GeV, a = 1/137.036 jlH], 
and the one- loop formula for a^J^^\ij) given in Eq. 0, with nj = 3 active quark flavors 
and Aqqj-, = 204 MeV |30|31j . Our default choice of renormalization scales is /i = m-p. 



where rriT = \Jpj- + 'im? is the transverse mass of the J/ip meson, and A = m. For 
the LO evaluation of the single- and double-resolved contributions, we use the LO set of 
photon PDFs from Gliick, Reya, and Schienbein (GRS) |30j, which are the only available 
ones that are implemented in the fixed-flavor-number scheme, with = 3, and iden- 
tify the factorization scale with /x. In want of NLO sets of J/V', Xcj-, and ip' MEs, we 
adopt the LO sets determined in Ref. [2Zj using the LO set of proton PDFs from Mar- 
tin, Roberts, Stirling, and Thorne (MRST98LO) [HI]. Specifically, (C"/'("^) fef^ 



' QXcQ 



3p(l) 



and 

were extracted from the measured partial decay widths of ip{nS) l^l^ 
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and Xc2 ^ 
and / 0^=0 



77 HE!, respectively, while (C"^("^) ^S^ ), (oV-CnS) s^j^ 



.(8)' 



Qil}{nS) 



^s[^^ ) were fitted to the transverse-momentum distributions of ip{nS) and Xcj 



inclusive hadropro duct ion ^5] and the cross-section ratio cr^^a/*^ 



Tevatron. The fit results for {^C"/'("^) 
so that the linear combination 



1^(8) 



and ( C"/'("^) 



(8) 



Xcl 



measured at the 



are strongly correlated, 



^il)(nS) 1^(8) 



3p(8) 



(19) 



with a suitable value of r, is quoted. Unfortunately, Eq. (fTHj) is sensitive to linear com- 



bination of (^V-CnS) 15-^8) \ /Q^inS) 



in Eq. ^T^i . In want of more specific in 



^Pg^^ ) that is different from the one appearing 
brmation, we thus make the democratic choice 

We now discuss the photon flux functions that enter our predictions for photoproduc- 
tion in the e~^e~ mode of TESLA. The energy spectrum of the bremsstrahlung photons 
is well described in the Weizsacker- Williams approximation (WWA) [33j by Eq. (27) of 
Ref. We assume that the scattered electrons and positrons will be antitagged, as was 
usually the case at LEP2, and take the maximum scattering angle to be 6'max = 25 mrad 
The energy spectrum of the beamstrahlung photons is approximately described by 
Eq. (2.14) of Ref. jHSl- R is controlled by the effective beamstrahlung parameter T, which 
is given by Eq. (2.10) of that reference. Inserting the relevant TESLA parameters for the 
y/S = 500 GeV basehne design specified in Table 1.3.1 of Ref. [33 in that formula, we 
obtain T = 0.053. We coherently superimpose the WWA and beamstrahlung spectra. 

At TESLA, prompt photons with scattering angles 7° < 6 < 173° are expected to 
be detectable, as can be read off from Fig. 2.4.2(a) in Ref. [23 ■ This corresponds to 



the pseudorapidity window \y^\ < 2.79. As a typical discrimination criterion against 
bremsstrahlung photons, we enforce the hardness condition pj. > Pr^mm with Pj'.min = 
2m = 3 GeV. This choice can be justified by studying the min variation (see Fig. [7| 
below). Unless otherwise stated, we always impose these acceptance cuts on pj. and y'^ in 
the following. 

With the NLO corrections to single- and double-resolved photoproduction yet to be 
evaluated, we are not in a position to present a complete phenomenological prediction 
that could be confronted with experimental data as it stands. Therefore, we refrain from 
presenting a full-fledged quantitative estimate of the theoretical uncertainties. However, 
a first indication of their size may be obtained by investigating the dependences of our 
numerical evaluation on the renormalization scales fi and A. For this purpose, it is suffi- 
cient to consider a typical kinematic situation. We thus choose as our reference quantity 
the differential cross section (Pa/dprdy at = 5 GeV and y = 0. 

In Fig. [3 the NLO result (solid line) is shown as a function of fi for mj'/4 < fi < ArriT, 
while A is kept fixed at its reference value. For comparison, the LO result (dashed line), 
which is of course /x independent, is also shown. As expected, we observe that the fi 
dependence of the NLO result is rather sizeable, reflecting the unscreened appearance of 
as{n) as an overall factor of the radiative correction. A partial compensation of this /x 
dependence can only occur at next-to-next-to-leading order, which is, however, beyond 
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the scope of this work. On the other hand, it is commonly beheved that the size of 
unknown higher-order corrections can be estimated by scale variations over judiciously 
chosen ranges of values. In this sense, we read off from Fig. El that the NLO result varies 
by ^32% over the interval mr/2 < fi < 2mT- We also observe that the QCD correction is 
negative for all values of /i considered here. This feature is scrutinized below, in connection 
with Figs. Hand El 

In Fig. ini the NLO result (solid line) is shown as a function of A for m/2 < A < 4m, 



while /i is kept fixed at its reference value. Since (C^ [^Sj^^J )^ (A) and (C^ [^P^^'l )^ (A) 
only enter Eq. (fTH|) at NLO, through dago and ddha, we keep them at A = m. Then, 
the A-dependent terms in Eq. (fT^ are just dao and 6op. In order to exhibit the partial 
compensation in A dependence between these two terms, we also include in Fig. |B1 the 
results that are obtained by only varying A in dao (dashed line) or 6op (dotted line) at 
a time. The dotted line is straight, reflecting the fact that 6op depends on A through a 
single logarithm. The theoretical uncertainty related to the A variation in the interval 
m/2 < A < 2m amounts to ^q\%- 

In the remainder of this section, we explore the phenomenological consequences of our 
analysis by studying the size and impact of the NLO corrections on the pj. > pj. ^^jj^, pt 
and y distributions of the cross section. In Fig. [71 we examine d'^a/dpx dy for px = 5 GeV, 
y = 0, and \y"'\ < 2.79 as a function of at LO (dashed line) and NLO (solid line). At 

LO, the prompt photon balances the transverse momentum of the J/ip meson, so that we 
integrate over a delta function located at pj^ = pt as long as pj^ < pt, which explains 
the step-function shape of the pj^ distribution. On the other hand, at NLO, there 
are events with an additional gluon-initiated hadron jet (j) to share the recoil transverse 
momentum of the J/tp with the prompt photon, so that, in particular, the latter attains 
access to the phase-space region where pj^ > pt- Of course, such events are suppressed 
by the smallness of as and, for pj, 3> pr, also by phase-space limitation. This explains 
the size and shape of the tail in the NLO result for pj. > pt- 

In Fig.jHJ we study d'^a/dpr dy (a) for ?/ = as a function of pr and (b) for px = 5 GeV 
as a function of y. In each case, the LO (dashed line) and NLO (solid line) results of 
direct photoproduction are shown. In Fig. jSl^a), the LO results of single-resolved (dotted 
line) and double-resolved (dot-dashed line) photoproduction are displayed as well. Notice 
that our analysis is only valid for finite values of pt] in the limit pt —>■ 0, additional 
IR singularities occur, which require a more sophisticated scheme of phase space slicing. 
Therefore, we do not consider px values below 2 GeV in Fig. IHl^a). 

As explained in the context of Fig. [3 we have px = Pt at LO, which explains the sharp 
threshold at pt = Pt in the dashed, dotted, and dot-dashed curves in Fig. jSI^a). By 
contrast, the corresponding threshold in the solid curve in Fig.jHl^a) is washed out by J/'ip'jj 
events. Comparing the LO results in Fig.jHfa), we observe that the direct-photoproduction 
contribution is overwhelming. It exceeds the single-resolved one, which is exclusively 
generated by the CO partonic subprocess 'y + g^cc ^S[^^ +7, by more than three orders 



3p(8) 



^ 3x10-3. 

For Pt ^ 6 GeV, the double-resolved contribution is chiefly produced by the CS partonic 



of magnitude, essentially reflecting the ME ratio (O-^/'^ ^S[^^ ) I lO^I'^ '^S\ 
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subprocess g + g cc ^Si^^ + 7, but it is nevertheless substantially suppressed by two 
factors of the gluon PDF of the photon and thus significantly undershoots the direct one, 
by approximately three orders of magnitude. At large values of px, the CO partonic 
subprocess q + q ^ cc ^Si'^ + 7 takes over the lead. This is typical for fragmentation- 
prone partonic subprocesses [T7j, which contain a gluon with small virtuality, = 4m^, 
that splits into a cc pair in the Fock state n = ^Si'^ and thus generally generate dominant 
contributions at pt ^ 2m due to the presence of a large gluon propagator. The relative 
importance of the three photoproduction modes is likely to be subject to change, at 
least at large values of pt, once the as-yet unknown NLO corrections for the single- 
and double-resolved contributions will be included. In particular, the NLO correction to 
the single-resolved contribution is expected to be sizeable because fragmentation-prone 
partonic subprocesses, which are absent at LO, start to contribute. In fact, a similar 
situation was encountered in Ref. [l^j . 

From Figs. IHl^a) and (b), we observe that the NLO result of direct photoproduction 
undershoots the LO one wherever the latter is non-zero, except at some edges of phase 
space. The magnitude of this reduction increases with the value of pt and decreases 
towards the forward and backward directions. This is nicely illustrated in Figs. IHl where 
the NLO to LO ratio K (solid line) is shown (a) for ?/ = as a function of pr and (b) for 
Pt = 5 GeV as a function of y, respectively. 

At first sight, the steady fall-off of the K factor as a function of pt in Fig. Efa) is 
surprising, and it is interesting to identify its origin. To this end, we break up the K 
factor into the contributions related to the cc Fock states n = ^S[^\ ^S^\^Plf^ as 



K=l + S -'SI 



3c(l) 



6 



10(8) 



+ 6 



3p(8) 



(20) 



where 6 



(1) 



10(8) 

On 



and 6 



3p(8) 



stem from the real 



10(8) 
^0 



and 6 



3p(8) 



is due to the virtual correction, and 6 

correction. Since the parts of the NLO correction that enter 6 
stem from different cc Fock states, they are separately gauge independent. However, the 
first and third of them are IR divergent, so that we are led to consider their MS-subtracted 
finite remainders. In order not to artificially introduce large logarithms, we choose the 
subtraction points to be yU = and A = m, respectively. In order to illustrate how the K 



and 1+6 



+6 



factor in Eq. ^lU^ is gradually built up, the results for 1-1-5 ^5*^^'' 
are also shown in Figs. IHfa) and (b), as the dashed and dotted curves, respectively. We 
observe from Figs. Efa) and (b) that the K factor is dominated by 6 '^Si'^ , while 6 



and 6 



3p(8) 



10(8) 

are of minor importance, which reflects the proportion of the respective 



MEs, (O^/^ ^Sl 



1^(8) 



we have (O^/'^ \^S[^^]) /M^i^ ^ 15 |27|. We also notice that 6 



3p(8) 



appearing as overall factors. In fact, 
is always positive. 



while the MS-subtracted quantity 6 



3p(8)- 



is negative for Pt^8 GeV. We conclude that 
6 ^5'}^'' is responsible for the striking fall-off of the K factor with increasing values of pt- 
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4 Conclusions 



Wc calculated the inclusive cross section of 77 J/ip+X^, where the system contains 
a prompt photon, in direct photoproduction for finite values of pr at NLO within NRQCD 
allowing for the J/ip meson to be promptly produced, and presented phenomenological 
predictions for TESLA in the e~^e~ mode of operation. At LO, 77 J/V' + pro- 
ceeds almost exclusively though direct photoproduction because this is a CS process. Its 
cross section is sizeable and its signal spectacular. These observations provided a solid 
motivation for our work. 

As for the real corrections, we employed the phase-space slicing method to demarcate 
the regions of phase space containing soft singularities from the hard regions, where the 
phase-space integrations were carried out numerically. We verified that the combined 
result is, to very good approximation, independent of the choice of the cut-off parameter 
6, over an extended range of values. We worked in dimensional regularization in con- 
nection with the MS renormalization and factorization schemes, so that our NLO result 
depends on the QCD and NRQCD renormalization scales fi and A, respectively. While 
the A dependence is formally cancelled up to terms beyond NLO, the /i dependence is 
unscreened and introduces an appreciable theoretical uncertainty. Our NLO result does 
not involve any factorization scale and is thus formally independent of single-resolved 
photoproduction and photon fragmentation. 

We found that the inclusion of the NLO correction leads to a substantial reduction in 
cross section, especially towards large values of px and in the central region of the y range. 
In fact, the K factor falls below 0.1 for px^ 12 GeV. These features could be traced to 
the (MS-subtracted) virtual correction, which makes up the bulk of the NLO correction. 

In order to complete the NLO treatment of 77 — > J/t/j+X^, we still need to evaluate the 
NLO corrections to single- and double-resolved photoproduction. Due to the occurrence 
of fragmentation-prone partonic subprocesses, these are expected to lead to a substantial 
enhancement in the case of single-resolved photoproduction. Nevertheless, it is unlikely 
that the dominance of direct photoproduction will be challenged at NLO. Then, also 
J/V' + X^ production in photoproduction at HERA and hadroproduction at the Tevatron 
can be described at NLO. This will help to provide a solid basis for an ultimate test of 
the NRQCD factorization framework. 
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Figure 1: Tree-level diagrams pertinent to the partonic subprocess 1^. 
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(a) 

Figure 2: One-loop (a) self-energy, (b) triangle, (c) box, and (d) pentagon diagrams 
pertinent to the partonic subprocess 
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(b) 

Fig. 121 (continued). 
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Figure 4: Tree-level diagrams pertinent to the partonic subprocess 
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Figure 5: Differential cross section (Pa/dprdy in pb/GeV of e+e~ — > e+e~J/'0 + Xj 
in direct photoproduction at TESLA with i/i = 500 GeV for a prompt J/t/j meson with 
Pt = 5 GeV and y = accompanied by a prompt photon with pj<>3 GeV and < 2.79. 
The LO (dashed hne) and NLO (sohd hne) results are shown as functions of /x. 
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Figure 6: Differential cross section dP'a / dpT dy in pb/GeV of e+e" e"'"e~J/'0 + 
in direct photoproduction at TESLA with = 500 GeV for a prompt J/^l) meson with 
Pt — 5 GeV and y — accompanied by a prompt photon withpj- > 3 GeV and < 2.79. 
The NLO result is shown as a function of A, (i) keeping A in the partonic cross sections 
fixed (dashed line), (ii) keeping A in (^O^ ^S'^^'' ^ (A) fixed (dotted line), and (iii) varying 
all occurrences of A simultaneously (solid line). 
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Figure 7: Differential cross section (Pa/dprdy in pb/GeV of e'^e~ — > e^e~J/'^ -\- 
in direct photoproduction at TESLA with ^/s — 500 GeV for a prompt J/'4> meson with 
Pt = 5 GeV and y = accompanied by a prompt photon with pj, > p^^^m ^^'^ Iv'^l < 2-79 
as a function of p^, j^^j^^ at LO (dashed line) and NLO (solid line). 
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Figure 8: Differential cross section (fa/dpT dy in pb/GeV of e+e" — >• e'^e~J/ip + at 
TESLA with ^/s — 500 GeV for a prompt J/ip meson accompanied by a prompt photon 
with > 3 GeV and \y'^\ < 2.79 (a) for y = as a function of pr and (b) for pt — 5 GeV 
as a function y at LO (dashed hne) and NLO (sohd hue) . In part (a) , the LO results of 
single-resolved (dotted line) and double-resolved (dot-dashed line) photoproduction are 
also shown. 
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Fig. IHl (continued). 
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Figure 9: QCD correction factor K of the differential cross section (Pa/dprdy in pb/GeV 
of e+e~ — > e+e~J/'0 + at TESLA with a/s — 500 GeV for a prompt J/ijj meson 
accompanied by a prompt photon with pj. > 3 GeV and \y'^\ < 2.79 (a) for y = as a 
function of pr and (b) for p^ = 5 GeV as a function y. The incremental shifts 1 + 6 ^S[^^ 
(dashed hues) and 1 + 5 \^S{^^] + 6 \^S^^^] (dotted hues) are shown separately. 
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Fig. ini (continued). 
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